Exploring Lattice Methods for Cold Fermionic Atoms 

Matthew Wingate a 

institute for Nuclear Theory, University of Washington, Seattle, WA 98195-1550 

There has been a surge of experimental effort recently in cooling trapped fermionic atoms to quantum degener- 
acy. By varying an external magnetic field, interactions between atoms can be made arbitrarily strong. When the 
S wave scattering length becomes comparable to and larger than the interparticle spacing, standard mean field 
analysis breaks down. In this case the system exhibits a type of universality, and J-W. Chen and D.B. Kaplan 
[1] recently showed how this system can be studied from first principles using lattice field theory. This poster 
presents the first results of exploratory simulations. The existence of a continuum limit is checked and the pairing 
condensate is studied as a function of the external source strength over a range of temperatures. Preliminary 
results show simulations can locate the critical temperature. 



1. INTRODUCTION 

At T = the physics of a dilute gas is charac- 
terized by just two length scales: the interparticle 
spacing n -1 / 3 = (-J^) 1 ^ 3 and the S wave scatter- 
ing length a. The range of interactions is very 
small compared to the average interparticle spac- 
ing R <C n -1 / 3 , so only 2-body S wave scattering 
is important. By varying an external magnetic 
field experimentalists can get a trapped Fermi gas 
with any scattering length, (e.g. [2]). 

Depending on the scattering length, the ground 
state of cold Fermi gases can be very different. 
In two extreme cases, theoretical expansions in 
(aln 1 / 3 <C 1 are valid. If there is a deep 2-fermion 
bound state (i.e. if < an 1 / 3 <C 1) then the 
fermion molecules form a Bose Einstein Conden- 
sate (BEC) at low temperatures. If there is no 2- 
body bound state (—1 -C art 1 / 3 < 0) then delocal- 
ized Cooper pairs condense in a Bardeen-Cooper- 
Schrieffer (BCS) supernuid. What is the ground 
state for fermions when \a\ is large? No perturba- 
tive expansion converges for lain 1 / 3 > 1. Ref. [1] 
is the only first principles approach for theoreti- 
cal calculation with lain 1 / 3 > 1. The case where 
| a | — > co is especially interesting since systems 
should then exhibit universal behavior. 

2. LATTICE THEORY 

Below we summarize the proposal of Chen and 
Kaplan [1]. Consider 2 species of nonrelativistic 



fermions, ip = 




, with a local 4-fermion 



interaction with finite chemical potential fi 
C = ^(^-e^-eo) - ^^(VV), 

The 2 species can be envisioned as spin degrees- 
of-freedom, although in the atomic gases they are 
two hyperfine states. 

In order to perform simulations we must make 
the Lagrangian quadratic in the fermion field by 
introducing an auxiliary real scalar field (f> so that 
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where we identify to 2 with — Cq 1 . 

Now one introduces a source for the fermion 
pairing, a necessary step to study spontaneous 
symmetry breaking in a finite volume 

C[J] = C + i(j^ T (T2V> + J*?CT2? T ) • (3) 

Combining tp and ijj into a 4-component Majo- 
rana fermion, ^f T = (ip T , ip(ia2) T ), the action can 
be written in terms of AVN t x AVN t matrices as 
^A^, where 

£)C? £)■ (4) 

K = K[4>] represents the fermion kinetic en- 
ergy and interaction terms, and it is the only 



2 



2VN t x 2VN t block to contain nontrivial space- 
time terms. The only coupling between species 
is through the itx 2 blocks. The partition function 
can be seen to have a real, nonnegative integrand: 



Z[J] = y"[d0][d*]exp(~ 
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(5) 



where K = K®h defines the VN t x ViV t matrix 
K . We can perform Monte Carlo simulations of 
(5) using Hybrid Monte Carlo. 

3. PARAMETER SPACE 

We have 3 lattice parameters, m 2 , /x, and M 
which control the physical quantities a, n, and 
T (ignoring finite V and setting N t fixed). 1/a 
is matched to the ratio m 2 /M at T, \i — by re- 
quiring 2-body scattering amplitude to match the 
leading term in the effective range expansion as 
in pionlcss nuclear effective field theory [3]. n is 
mainly controlled by tuning fi. In a nonrelativis- 
tic theory, the fermion mass can be scaled away by 
rescaling the time coordinate, so it can be inter- 
preted as the anisotropy. For fixed N t , increasing 
M shortens the temporal length, increasing T. 

4. CONTINUUM LIMIT 

Some parts of parameter space do not describe 
the physics of a dilute Fermi gas. In a continuum 
theory with a local attractive interaction, the vac- 
uum is unstable. Since we are not interested in 
physics at short distances, we simply have to reg- 
ulate the theory. In our case, the lattice automat- 
ically provides a cutoff (see Fig. 1). 

The continuum limit is defined as the limit 
where physical correlation lengths a,n -1 / 3 , be- 
come infinitely large compared to the lattice spac- 
ing b. This limit can only be taken using simula- 
tions where one finds n6 3 < 1 at the same time 
that 1/a is tuned to 0. 

Figure 2 shows n <C 1 (in lattice units) at 
1/a = (and ji = 0) which is necessary and 
sufficient for a continuum limit to exist, n is 
nonzero there due to finite V,N t and J effects. 
As 1/a — > 1, a rise in n is observed. It would 
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Figure 1. Mean field illustration of the instability 
at strong coupling. The solid line corresponds to 
a coupling which is small enough in lattice units 
for the ground state to be the physical one, the 
dashed line corresponds to a coupling so strong 
that the ground state favors a fermion at every 
site even at fi = 0. 



be satisfying to be able to increase 1/a enough so 
that we observe the artificial first-order transition 
to the filled-lattice state, but the HMC algorithm 
is becoming less efficient in this direction. For 
larger values of 1/a than plotted, 6000 HMC tra- 
jectories were not enough to equilibrate the lat- 
tice. The breakdown of the HMC algorithm is 
known from simulations of the Hubbard model. 
Indeed the physics in the uninteresting (from our 
perspective) region of parameter space is exactly 
that of the 3D attractive Hubbard model. Having 
established the existence of a continuum limit in 
the interesting part of parameter space, searching 
for the uninteresting region is not a priority. 

5. CRITICAL TEMPERATURE 

At low temperatures we expect condensation of 
fermion pairs, breaking the global U(l) symme- 
try of fermion number conservation. The order 
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Figure 2. Number density vs. inverse scattering 
length (both in lattice units) computed on 6 3 x 12 
lattices with J = 0.05 and a = 0. Since n <C 1 in 
lattice units, we can conclude the ground state is 
the physical one and the continuum limit exists. 



parameter is 

s = - ^(V> T <7 2 V + V^2? T ) . (6) 

Figure 3 shows S as a function of J for several 
values of M. Extrapolating to J = linearly 
indicates the vanishing of the order parameter as 
M increases. The validity of the linear extrapola- 
tions is under study, since we expect finite volume 
errors to increase as J decreases [4]. 

6. CONCLUSIONS 

The results of this exploratory study show this 
lattice theory can describe the physics of dilute 
Fermi gases. These simulations show a contin- 
uum limit exists: it is defined by bn 1 ^ 3 ,b/a — ► 
with fixed an 1 / 3 . Study of the critical tempera- 
ture has begun. The first point on the critical line 
has been located in the space of lattice parameters 
(m 2 , M). The critical line can then be mapped to 
the physical space (1/a, T). For precision results, 
one will have to repeat these steps for increasing 



Figure 3. Order parameter vs. pairing source J 
(both in lattice units) computed on 8 3 x 16 lat- 
tices with /i = 0.4 and m 2 = 0.1456. Increasing 
M corresponds to increasing T; linearly extrap- 
olating the order parameter to J = indicates 
U(l) restoration close to M = 1.4. 



volumes and decreasing lattice spacings. How- 
ever, even with finite size and spacing effects, un- 
derstanding the behavior of T c as a function of 
an 1 / 3 will be an important accomplishment. 
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